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We propose a method for measuring entangled vibronic 
quantum states of a trapped atom. It is based on the nonlin- 
ear dynamics of the system that appears by resonantly driving 
a weak electronic transition. The proposed technique allows 
the direct sampling of a Wigner-function matrix, displaying 
all knowable information on the quantum correlations of the 
motional and electronic degrees of freedom of the atom. It 
opens novel possibilities for testing fundamental predictions 
of the quantum theory concerning interaction phenomena. 

PACS numbers: 03.65. Bz, 42.50.Vk, 32.80.Pj 



I. INTRODUCTION 

Entanglement is one of the most striking aspects of 
quantum mechanics In classical physics, two inter- 

acting systems retain their individuality during the in- 
teraction process and become completely independent of 
each other after their coupling has been switched off. By 
way of contrast, quantum theory predicts a completely 
different behaviour of interacting systems. When two 
quantum systems are brought to interact, their identities 
become in the course of time more and more entangled, 
so that a state-vector description of each system is in 
general precluded. They build an entangled composite 
system, whose state-vector cannot be separated into a 
product of the states of the subsystems. More surpris- 
ing, the entanglement is preserved even after switching 
off the interaction, so that a measurement on one system 
will affect the other one. The individual systems can be 
characterised by (reduced) density matrices. However, 
even when the reduced density matrices of both subsys- 
tems are known, important information on the physics of 
their interaction and the related entanglement is lost. 

Recently, entangled quantum states of an electronic 
system and a harmonic oscillator have been realized with 
trapped atoms |^ and in cavity QED Q . The signatures 
of the entanglement have been partially demonstrated by 
quantum measurements. The determination of the full 
quantum statistical information of such quantum states, 
however, requires new types of quantum measurements. 
Moreover, it is interesting in this context that with both 
types of systems the Jaynes-Cummings model |5] and 
its nonlinear trapped-atom counterpart could be real- 
ized 1^,^ . In the Jaynes-Cummings interaction the quan- 
tum correlation between the interacting systems has been 
predicted to show surprising features |y. In all these 
cases methods for determining the full information on 
the quantum state of entangled systems are desired for 



getting more insight into interaction phenomena in the 
quantum world. 

Several approaches have been proposed for determin- 
ing the complete information on the quantum state of 
single quantum systems. In particular, homodyne to- 
mography has been established for optical fields pO[ . 
Further simplifications of the method have been intro- 
duced, including the direct statistical sampling of the 
density matrix in field-strength (quadrature) represen- 
tation pT|-p^, photon- number representation [Q, and 
symplectic tomography |l5[| . Tomographic methods have 
also been realized for molecular vibrations and pro- 
posed for reconstructing the motional quantum state of 
trapped atoms Alternatively, phase-space distribu- 
tions can be determined by measuring the number statis- 
tics of the quantum state of interest, after introducing ap- 
propriate coherent displacements |]l8|. A method of the 
latter type has recently been used to reconstruct the mo- 
tional state of a trapped atom [l9| ] . For quantum systems 
undergoing a Jaynes-Cummings type dynamics H, such 
as high-Q cavity fields |^,|| or trapped atoms pi|,p|,[7t , 
various methods have been considered to study the quan- 
tum statistics of the bosonic subsystem (the field or the 
center-of-mass motion), cf. ||2^- |2^Jl9| . All these meth- 
ods, however, give no insight into the entanglement of 
quantum systems. 

In this contribution we propose a method for determin- 
ing entangled vibronic quantum states of a trapped atom. 
In Sec. H of this paper we introduce a Wigner-function 
matrix, which contains the full information on the com- 
posite system under study. After a brief discussion of its 
properties we show that it can be determined via mea- 
surements of entangled, motional number statistics. In 
the following section we present a concrete scheme for 
measuring those quantities. It relies on the nonlincari- 
ties appearing in the motional dynamics of the trapped 
atom interacting with laser light, which have been pre- 
dicted for atomic localisations beyond the Lamb-Dicke 
regime and confirmed in recent experiments [Q. By 
monitoring the electronic dynamics, combined with co- 
herent displacements of the motional subsystem, one can 
perform a direct sampling of the Wigner-function matrix 
of the vibronic state. In Sec. |^ wc present a numerical 
simulation of the reconstruction of the Wigner-function 
matrix and briefiy discuss the practical aspects of our 
method. A summary and some conclusions are given in 
Sec.0. 
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II. CHARACTERISATION OF ENTANGLED 
QUANTUM STATES 

The motional quantum state of a trapped atom is usu- 
ally described by a density operator p, which is an oper- 
ator in the Hilbert space Ticm of the motion of the atom 
in the trap. To get a phase-space description of the mo- 
tional state of the atom one is lead to quasiprobability 
distributions such as for example the Wigner function. 
This function can be obtained from the density operator 
as follows [Esl^Pl 



W{a) = ((5(a-a)) = Tr p5{a-a) 



(1) 



where the operator- valued delta-function 5{a — d) is the 
Fourier-transform of the displacement operator D{a) = 
eyiY>{aa' — a* a) and reads as pSf 



5{a-a) = -y d2^I)(e)e"«*-"*« 

IT 



(2) 



Equation (Q) can be inverted to get the density operator 
out of the Wigner function 



(3) 



p= d^aW{a)5{a-a), 



so that the knowledge of W{a) is equivalent to the knowl- 
edge of the motional density operator p. 



A. Wigner-function matrix 

Even though the Wigner-function 1/1^(q;) displays all 
the obtainable knowledge about the motional state of the 
atom, it gives no information on its electronic degrees of 
freedom. For a complete description of the atomic state 
— i.e. including the electronic degrees of freedom — one 
has to generalise the concepts leading to the definition of 
the Wigner-function. For this purpose we introduce the 
density operator g of the whole system, which is now an 
operator in the product-space of motional and electronic 
Hilbert spaces TYcm ® 'Hd- The reduced density opera- 
tors for the motional and electronic subsystems can be 
obtained by taking the trace over the Hilbert spaces of 
the electronic system (el) and the center-of-mass motion 
(cm), respectively, 



Troi [q], ct = Trcm [q\ ■ 



(4) 



We now define a Wigner-function matrix describing the 
complete quantum state of the trapped atom by general- 
ising Eq. (0), 



Wij{a) = {S,j{a-a)) = Tr gSij{a-a) 



(5) 



where Sij{a—a) is now an operator in the product-space 



5ij{a — a) = Aji S{a — a). 



(6) 



Here we have used the electronic flip-operators Aji = 
\j){i\, which describe transitions from the electronic state 
|z) to the state In close analogy with the case of the 
motional sub-system (||) , the density operator of the com- 
posite system can be obtained from this Wigner-function 
matrix as 



g — J (fa Wij (a) Sji{a — a), 



(7) 



showing that the complete information on the quantum 
state of the atom is contained in the Wigner-function 
matrix Wij (a). 



B. Properties of the Wigner-function matrix 

Let us now discuss some properties of the Wigner- 
function matrix Wij (a). It is easy to see that its di- 
agonal elements Waia) represent the Wigner functions 
of the conditioned density operators ga correlated to the 
electronic state |z). More precisely, this means that if one 
would measure the electronic state of the trapped atom to 
be in state \i), the corresponding (unnormalised) condi- 
tioned Wigner function of the motional sub-system would 
be Wii{a). The norm of this conditioned Wigner function 
is simply the occupation probability an of the electronic 
state \i), whereas the norm of the off-diagonal elements 
Wij{a) (i^j) represents the electronic coherence 



(fa Wij (a) = Tr 



A-ji g 



(8) 



Moreover, from Eqs. (||) and (^ it is seen that the 
Wigner-function matrix is hcrmitian 



W,jia) = W*,ia), 



(9) 



so that in order to have the full information on the quan- 
tum state under consideration, it is sufficient to know 
the real- valued diagonal elements of the Wigner-function 
matrix Wii{a) and real and imaginary part of the off- 
diagonal elements Wij{a) {i<j). 

If one measures the electronic state to be in an arbi- 
trary superposition 



(10) 



the conditioned density matrix of the motional sub- 
system is 



P 



'j Qij- 



(11) 



He 



) Hei and is defined by 



The Wigner function of this motional quantum state can 
now be represented in terms of the Wigner-function ma- 
trix (i) 



2 



(12) 



This result shows that the ofF-diagonal elements of the 
Wigner-function matrix, which are in general complex- 
valued, contain the information on the electronic coher- 
ences and the corresponding motional states, which are 
entangled with these coherences. 

On the other hand, if one has no information about 
the electronic state of the atom, the density operator of 
the motional sub-system is given by 



Tr, 



(13) 



The Wigner function of this reduced density operator 
reads as 



M/(a) =Ew,,(a), 



(14) 



and is simply the trace of the Wigner-function matrix. 

In Figure 1 we show an example for the Wigner- 
function matrix of an entangled state corresponding to 
the so called "Schrodinger-cat" states. The quantum 
state considered here has been recently prepared experi- 
mentally with a trapped atom by Monroe et al. ||^. It is a 
quantum superposition of two motional coherent states, 
I ±/3), of amplitudes /3 and — /3, entangled with the upper 
and lower electronic states, respectively. This state can 
be given by the following expression 



i( 1/3)12) 



x/2 



The Wigner-function matrix for this entangled quantum 
state can be explicitely given by 



W22(a) 



TT 

2 



■ exp(— 2|a 

■ cxp(— 2|a 



(16) 



exp(-2|a|^)exp[2iIm(a/3*) 



It is clearly seen, that Wii(q;) and W22(a) represent the 
coherent states with amplitudes —(3 and (3, respectively, 
whereas Wi2{oi) contains the information on the elec- 
tronic coherence and the quantum interference effects in- 
herent in the entangled state under study, see Fig. 1. 

If the quantum state under consideration can be writ- 
ten as a product state 



Q — p® a, 



(17) 



the corresponding Wigner-function matrix reads as 

W,,{a)^a,^W{a). (18) 

States of this type contain no entanglement between the 
motional and electronic degrees of freedom. The corre- 
sponding Wigner-function matrix Wij {a) is of the same 
shape for all indices determined by the motional 
Wigner function iy(a) and weighted by the electronic 
density matrix elements dij . 





Re{a) 



Im(a) 



(15) Re[Wu{a)] 




Re{a) 



Im(a) 



Im [Wuia)] 



0.2 




Im(a) 



FIG. 1. Wigner-function matrix Wjj (a) of the entangled 
vibronic quantum state \ip) [cf. Eq. (M)] with /3 = 2. 
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C. Representation by displaced number statistics 

For determining the Wigner-function matrix Wij (a) 
from measured data we note that it can be represented 
in terms of matrix elements of the coherently displaced 
density operator g, where only the motional diagonal ele- 
ments are needed. To show this we make use of Eqs. §), 
together with Eq. (||). By taking the trace over the 
electronic sub-system, the Wigner-function matrix can 
be written as 



W,j{a) = - Tr [ft, D{a) {-lf^D^{a) 



(19) 



where gij = {i\g\j) is still an operator acting on Ticm- By 
using the cyclic property of the trace and performing the 
trace in terms of number-states of the harmonic vibration 
in the trap, the Wigner-function matrix reads as 




|2) 



W,,{a)^-J2{-irgl^a). 



(20) 



n=0 



Here we have made use of the coherently displaced den- 
sity operator 



g{a) ^D\a)gD{a). 



(21) 



For representing the Wigner-function matrix we only 
need the diagonal elements of g{a) with respect to the 
motional degree of freedom 



(22) 



For convenience we will call these quantities displaced, 
entangled number statistics. 

To get the full information on the entangled quantum 
state, it is sufficient to measure the displaced, entangled 
number statistics (P2|), i n generalisation of the method 
presented in Ref. [[18| , |19[ . The initial coherent displace- 
ment can be realized by applying a radio-frequency field, 
as has been done in the experimental determination of 
the quantum state of the motional subsystem ||l9| . In the 
following we present a measurement scheme for determin- 
ing the displaced, entangled number statistics g^p{a). 



III. MEASUREMENT SCHEME 

The basic scheme consists in a three-level electronic 
system of V-type which is superimposed by the energy 
levels of the motion of the atom in the harmonic trap po- 
tential, see Fig. 2. A weak electronic transition |1) ^ |2) 
is the one of interest with respect to its entanglement 
with the motional degree of freedom. The dynamics of 
this transition is monitored with very high quantum effi- 
ciency by testing a strong, auxiliary transition 11) ^ |3) 
for the appearance of resonance fluorescence |26| . We as- 
sume that the weak transition is driven in the resolved 
sideband limit. When the laser is tuned on resonance 
with the vibrationless line, the corresponding interaction 
Hamiltonian (in the interaction picture) reads as p,CT 



FIG. 2. Three-level electronic system of the trapped atom 
for the measurement of the displaced, entangled number 
statistics. The weak electronic transition |1) <-> |2) is tested 
by probing the strong transition |1) <-> |3) for resonance fluo- 
rescence. 



iJint = -hnf{a''a)Ai2 + h.c. 



(23) 



Here f2= exp{i(p) is the Rabi frequency of the laser in- 
teracting with the weak electronic transition and ip is the 
phase of the laser field |^^. The operators a and are 
the annihilation and creation operators of a vibrational 
quantum, respectively. The function /(a^a), describing 
the nonlinearities in the vibronic coupling, is given in its 
normally ordered form by 



E 

k=0 



k„2k 



(kiy 



(24) 



with r] being the Lamb-Dicke parameter, characterising 
the localisation of the atom in the trap. The interac- 
tion Hamiltonian ( p3| ) fulfils the condition [n, H-mt] =0 of 
back-action evasion for the vibrational number operator 
n = a^a, so that the latter is a constant of motion. Con- 
sequently, this interaction couples only between vibronic 
density matrix elements g"^^ having the same motional 
indices {m,n). This renders it possible to determine the 
entangled motional number statistics g™ of the initial vi- 
bronic state by monitoring the dynamics of the electronic 
transition |1) ^ |2). 

One possibility to obtain the matrix elements g™{a) 
is the measurement of the time dependence of the oc- 
cupation of state 1 2). This dynamics is completely de- 
termined by those matrix elements. However, their re- 
construction requires some effort of data analysing by 
Fourier or other appropriate techniques, as considered 
for the Jaynes-Cummings dynamics |22 2^ , p^ ,p[| . 

In the following we will deal with an alternative 
method that directly yields the quantities desired in 
Eq. (|2^) . It is related to the scheme proposed in Ref. |2^] 
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for the quantum nondemolition measurement of the mo- 
tional energy of a trapped atom. After a well controlled 
interaction time of the laser resonant to the weak elec- 
tronic transition, a laser pulse is applied on the strong 
transition in order to probe for the appearance of reso- 
nance fluorescence. The occurrence of fluorescence de- 
tects the atom in the electronic ground state |1) and 
its absence in the excited electronic state |2). We will 
focus our attention to the no-fluorescence events, since 
they preclude any disturbance of the motional state via 
light scattering by the strong transition. Alternating se- 
quences of light pulses on the weak and the strong tran- 
sition allow to reduce the motional state to a Fock state, 
which is predetermined by time control of the pulse se- 
quence on the weak transition. Eventually, the desired 
information on the entangled number statistics g^"{a) 
is directly given by the probability of realizing such a 
sequence of interaction-probe cycles, for each predeter- 
mined Fock state \n). 

Assume that, subsequent to the coherent displacement 
of the motional subsystem, the atom is probed for flu- 
orescence on the strong transition after an interaction 
time T of the laser with the weak transition. Provided 
the atom is detected in the electronic state |2) (no fluo- 
rescence), the density operator of the system reduces to 
the (unnormalised) operator 

e(-^)(r) = |2)(2|®p(-<i)(T), (25) 

where /5*^"'*^(r) = (2|^(t)|2) is the corresponding (unnor- 
malised) density operator of the motional subsystem re- 
lated to the electronic state |2). In view of the interac- 
tion Hamiltonian (p^), the diagonal elements of p*^™'*)(t) 
in number-state representation are given by 

P^nn'Kr) = ,^,"(a)cos2 (^^) + e™(a)sin2 (^^^^ 

+ Im[f?™'(a)e-^'^]sin(a.T), (26) 

where r2„ = |i7|L„ (77^) exp(— 77^/2) are the nonlinear vi- 
bronic Rabi frequencies, L„{x) being Laguerre polyno- 
mials. 

Let us consider the effect of k of these interaction- 
probe cycles with interaction times ri, . . . , r^, each one 
accompanied by no fluorescence |^^. The resulting (un- 
normalised) motional number statistics is conditioned on 
the times ti , . . . , , at which the reductions to the state 
1 2) occurred, and reads as 

pt'Kr,, . . . ,n) = n cos^ (^) Pi'^Hn). (27) 

The probability P{Tk, . . . ,ti) to obtain such a sequence 
of interaction-probe cycles is given by the trace of 

P(rfc,...,ri) = ^p(-d)(r,,...,ri), (28) 

n=0 



and can be experimentally determined by repeating this 
procedure many times and counting the number of times 
one was successful in obtaining such a sequence. 

An adequate choice of the interaction times allows 
to map pnrf'\Ti) onto the probability P(rfc, . . . , ri). For 
mapping the particular element p^m{Ti), it is useful to 
choose the interaction times r2, . . . , as 

T2,...,Tk = —p, p=l,2,..., (29) 

where p can be set to a different integer at each inter- 
action cycle. From Eq. ( p7| ) it is seen that after an ap- 
propriate number k > fcmin of interaction-probe cycles 

Pn!f \Tk, • • • , Ti) reduces to 

p^Z''Hn,...,n)=6nrapt''rnHT^). (30) 

The other matrix elements {nj^ m) are suppressed by 
the product of cosines in Eq. (|2^) , provided that differ- 
ent motional number states are efficiently discriminated 
by the vibronic Rabi frequencies ri„ |2^. It is easy to 
see [cf. Eqs (11), (1^)] that in this case the probabil- 
ity P(rfc, . . . , Ti) reflects the quantity p^miri) and that 
further interaction-probe cycles will not change anything. 
Since the interaction times are specified from the be- 
ginning, the minimum number fcmin of cycles needed to 
complete the mapping process can be evaluated. 

The number statistics p^mm (ti ) obtained in this man- 
ner can be used to determine the entangled, motional 
number statistics g^"^{a) of the displaced initial vibronic 
state via Eq. ([2^). For this purpose, three different 
choices of the first interaction time ri are desired: 

(i) Ti = 0: After detecting the atom in the excited 
electronic state |2) (no detection of fluorescence), 

Prnm 

(ti) is given by 

ptt\ri)^Q'^2^{c^)- (31) 

(ii) Ti —Ti/VLm'- For this choice the diagonal elements 
of the motional density matrix are given by 

pt:^\T,)^Q\T{a)- (32) 

(iii) Ti =7r/(2ri„j): In this case one arrives at 

+Im [QXT{a)e~^^] . (33) 

Choosing two laser phases, (^ = and 7r/2, one 
can obtain the imaginary and real part of Qi2^{oi), 
respectively, by subtracting one half of the out- 
comes from (i) and (ii). 

Eventually, any positive event (series of no fluorescence) 
recorded in this manner can be stored in the computer for 
a direct sampling of the Wigner-function matrix (to be 
normalised by the number of trials) according to Eq. ([2^). 
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IV. NUMERICAL SIMULATIONS 



Wu(a 



0.3 



. 2 



0.1 



W-22(a 



A typical example for an entangled vibronic quantum 
state of a trapped atom is the state given in Eq. (15). 
In Fig. 3 we show a simulation of all steps of the scheme 
for the determination of the Wigner-function matrix of 
the state \tp) according to the method proposed above. 
The maximum pulse length of the laser resonant to the 
weak transition is of the order of 6/iS for a typical Rabi 
frequency of ri/27r = 500 kHz [Q. For interaction times 
of this order the electronic relaxation of the weak tran- 
sition is negligible. The time needed for the data ac- 
quisition of a point in phase space is approximately 30 
seconds, so that the complete sampling of one matrix el- 
ement of the Wigner-function matrix Wij{a) as shown 
in Fig. 3 would take about 3 hours. This is comparable 
with the data-acquisition time needed for the reconstruc- 
tion of the quantum state of the motional subsystem, by 
measurement of the time evolution of the electronic inver- 
sion jl9|. The advantage of our scheme, however, consists 
in the fact that the desired number statistics is directly 
observed, with no need for inverting systems of linear 
equations with noisy data input. 

The diagonal elements ^22(0) and 1^11(0!) simply rep- 
resent the Wigner functions of the coherent states \f3) and 
| — /3) respectively. The off-diagonal elements Wij{a) con- 
tain the more interesting information on quantum inter- 
ferences inherent in the entangled state (Jig) . Comparing 



with the theoretical result of Wij{a) given in Fig. 1, a Re [VFi2(ct)] 



good agreement between theory and simulation is found. 
Note that the entanglement is reflected by the fact that 
the Wigner-function matrix is different for different elec- 
tronic indices. If one would measure the Wigner function 
of the motional sub-system, iy(a) — J^i W^m(Q!)i the re- 
sult would be nothing but an incoherent superposition of 
two coherent states, completely concealing the nonclassi- 
cal character of the atomic state. That is, in the chosen 
example the nonclassical nature of the system is mani- 
fested by the entanglement. Our measurement scheme is 
directly suited to demonstrate these features. 





Re(a) 



Im(a) 




Re(Q!) 



Im(a) 




Im(a;) 



V. SUMMARY AND CONCLUSION 

In conclusion we have proposed a measurement tech- 
nique for obtaining the full information on entangled vi- 
bronic states of a trapped atom. First of all, this re- 
quires a characterisation of the full quantum state of the 
composed system which can be related to an appropri- 
ate measurement scheme. For this purpose we define 
a Wigner-function matrix which has the character of a 
density matrix with respect to the electronic degrees of 
freedom, and of a Wigner function for the motional sub- 
system. This Wigner-function matrix is readily related 
to the displaced, entangled motional number-statistics of 
the atom. 

For the determination of the Wigner-function matrix 
we have presented a method that allows the measurement 
of the displaced, entangled motional number-statistics. It 




FIG. 3. Wigner-function matrix Wij(a) of the entangled 
vibronic quantum state [Eq. (p^)] with /3 = 2. Each matrix 
element (0; a) has been numerically sampled with 1000 
trials and a sequence of fc = 30 interaction-probe cycles on a 
25 X 15-grid. The Lamb-Dicke parameter is ?7 = 0.1. 
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is based on the nonlinearities appearing in the vibronic 
couphng of a resonantly driven, weak electronic transi- 
tion. This coupling fulfils the back-action evasion crite- 
rion for the number of motional quanta. It eventually 
allows, in combination with electronic-state reductions 
by fluorescence measurements, the direct statistical sam- 
pling of a Wigner-function matrix displaying the entan- 
glement between the motional and the electronic degrees 
of freedom of the atom. The method proposed here opens 
new perspectives for experimental demonstrations of fun- 
damental properties of interacting quantum systems. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 
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